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Abstract: In this paper, a new technique is used to an 

unsupervised learning image classification based on 
integration between neutrosophic sets and optimization linear 
programming. Neutrosophic sets are used to segment the 
image into three main components namely objects (0), edges 
(E), and Background (IS). The neutrosophic image 

components ( 0,E,B ) are corresponding to the 
neutrosophic sets components (T , /, F). The components of 

neutrosophic image valued in [ 0,1 + ] are representing the 

association intensities degree of pixel for each image 
components. Netrosophic image components are contributed 
to solving one of the important problems in image 
classification known as "overlapping" within cluster. While, 
the problem of overlapping between clusters is solved by using 
optimization linear programming. 

Key words: Neutrosophic set, image classification, linear 
programming optimization 

I. INTRODUCTION 

Since several decades, the world is witnessing a remarkable 
development in the science of computer vision. The 
principle of computer vision is based on deal with the 
images and methods of treatment. Hence the interest of 
researchers in the computer vision with image processing, 
which is concerned, in essence, on the methods and many 
different algorithms. Among these algorithms are image 
classification algorithms. 

Classification is the field devoted to the study of methods 
designed to categorize data into distinct classes. This 
categorization can be divided to distinct labeling of the data 
(supervised learning [1]), division of the data into classes 
(unsupervised learning [2]), selection of the most 
significant features of the data (feature selection [3]), or a 
combination of more than one of these tasks [4] . 
Unsupervised image classification (UIC) starts by 
partitioning the image data into groups (or clusters). The 
classes in UIC are unknown, according to similarity 
measure, groups of image data can be compared with 
reference to data by an analyst [5]. UIC can be categorized 
into two main groups namely Hierarchical [6] and 
Partitional [7] algorithms. 

In hierarchical clustering algorithms (HCA) a sequence of 
clustering with each clustering being a partition of the data 
set are showing as a tree [8]. HCA is characterized by two 
advantages, first the number of classes does not need be 



specified a priori and the others they are independent of the 
initial condition. However, HCA is suffers from be a static 
algorithm and its inability to solve the overlapping clusters 
problem [9]. HCA are divided according to the clusters 
construction methods or according to the similarity 
measure. For methods construct the clusters by recursively 
partitioning the instances in either a top-down or bottom-up 
fashion. These methods can be subdivided as 
agglomerative [10] and divisive [11] methods. Whereas, 
the merging or division of clusters is performed according 
to some similarity measure, chosen so as to optimize some 
criterion (such as a sum of squares). The hierarchical 
clustering methods could be further divided according to 
the manner that the similarity measure is calculated [12], 

On the other hand, partitional clustering algorithms (PCA) 
are based on image data set segmentation into a specified 
number of clusters. PCA can be treated as an optimization 
problem as a result of reliance on the square error function 
to minimize certain criteria. Both HCA and PCA 
algorithms are participate in advantages and drawbacks. 
There are two categories from PCA namely Iterative [13] 
and non-iterative [14] algorithms. K-means algorithm [15] 
is the most widely used in iterative partitional algorithms. 
The basic idea for k-means algorithm is to find a clustering 
structure that minimizes a certain error criterion which 
measures the distance of each instance to its representative 
value. The most well-known criterion is the Sum of 
Squared Error (SSE) [16], may be globally optimized by 
exhaustively enumerating all partitions, which is very time- 
consuming, or by giving an approximate solution using 
heuristics. Another partitioning algorithm, which attempts 
to minimize the SSE is the k'-medoids [17] or partition 
around medoids (PAM) [18]. 

Lillesand and Kiefer [19] presented a non-iterative 
approach to unsupervised clustering with a strong 
dependence on the image texture. Researches [20-21] have 
shown that the iterative algorithms are more efficient than 
its counterpart non-iterative, where it does not rely too 
much on data points order. 

There are other unsupervised classifications methods are 
used recently represented in Density-based Methods [22] 
which assume that the points that belong to each cluster are 
drawn from a specific probability distribution. Model- 
based Clustering Methods [23], these methods attempt to 
optimize the fit between the given data and some 
mathematical models. Unlike conventional clustering, 
which identifies groups of objects; model-based clustering 
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methods also find characteristic descriptions for each 
group, where each group represents a concept or class. The 
most frequently used induction methods are decision trees 
[24] and neural networks [25]. Grid-based Methods [26], 
these methods partition the space into a finite number of 
cells that form a grid structure on which all of the 
operations for clustering are performed. The main 
advantage of the approach is its fast processing time [27]. 
Soft-computing Methods, In addition to neural networks, 
there are some methods that belong to soft computing 
methods such as Fuzzy Clustering [28], Evolutionary 
Approaches for Clustering [29] and Simulated Annealing 
for Clustering [30]. 

In this paper, a new an unsupervised image classification 
technique is used based on neutrosophic sets [31] and 
optimization linear programming [32]. Neutrosophic set 
(N s S) considered a part from neutrosophy theory which 

interested by studies the origin, nature and scope of 
neutralities, as well as their interactions with different 
ideational spectra. The idea of neutrosophy theory depends 

on event or entity, where between an idea (Aj and its 
opposite ( Anti — A'j , there is a continuum power 

spectrum of neutralities ( Neut — Aj . Truth value (T), 

indeterminacy value (I) and falsehood value (F) are 
represented neutrosophic components referring to 
neutrosophy, neutrosophic logic, neutrosophic set, 
neutrosophic probability, neutrosophic statistics [33]. In 
neutrosophic set, the indeterminacy is quantified explicitly 
and the truth-membership, indeterminacy-membership and 
falsity-membership are independent. The neutrosophic set 
is a generalization of an intuitionistic set [34], classical set 
[35], fuzzy set [36], paraconsistent set [37], dialetheist set 
[38], paradoxist set [39], and tautological set [40]. 

Linear programming is constrained optimization, where the 
constraints and the objective function are all linear. It is 
called "programming" because the goal of the calculations 
help you choose a "program" of action [41]. The linear 
programming model, for neutrosophic image classification 
problem, involves on two main parts called constraints and 
objective function. Constraints are describing the query 
images as lower and upper weights for neutrosophic query 
image components. On neutrosophic image clustering 
classification to be maximized a linear objective function 
means that categorization of similar images in clusters with 
out overlapping within or between clusters. 

The rest of the paper is organized as follows: section 2 
presents general framework for proposed technique, 
neutrosophic image processing is given in section 3. feature 
extraction for neutrosophic image is presented in section 4. 
Neutrosophic image cluster is illustrates in section 5. 
Section 6 illustrates accuracy evaluation for the technique. 
Section 7 presents experimental results to illustrate the 
efficiency of the algorithm. Section 8 concludes the paper, 
and outlines future research. 



II. GENERAL FRAMEWORK 
This paper presents a novel system to image clustering 
namely Optimization neutrosophic image classification 
system (ONsICS). As shown in figure 1, ONsICS consists 
of two techniques are neutrosophic image processing and 
optimization image clustering. Neutrosophic image 
processing is used to convert gray image to enhanced 
binary image (EBI) based on object, edge and background 
of image components. Each image can be represented as 
neutrosophic components (T, I, F) and stored the extracted 
image components feature as a vector in database. All 
similar image features are gathered together in a one 
category by using neutrosophic image clustering (NsIC) 
technique. Image clusters are optimized by using linear 
programming to solve image overlapping problem as 
shown in figure 1 . 



IIL NEUTROSOPHIC IMAGE PROCESSING: 



Let I img be a Universe of discourse represents image and I comp is a 



set of I img represents image components (as object, edge, 
background) which is composed by bright pixels. Aim of 
the neutrosophic image domain (N S D) is transferring 

image I img to neutrosophic domain by describing the pixel 
by three membership sets T, I and F as 
P Ns (T,I , F ) [43]. The pixel can represents as: 

P Ns (e J ) = Me j\ 1 (e j\ Me ./')} where , 

T(i, jj is the probability belonging to white pixels set. It is 
defined as: 



T(i,j) = M 



(e j)~ 



g r 



§ max & min 

1 t+% j+% 

g( uj )= — X 

wxw ,n=i-y 2 n=j-y 2 

Where, g(i, j j is the local mean value of pixels of the 
window. 

7 (e j) is indeterminate set. It is defined as: 






S „ -S. 



max mm 



S(i, j ) = ahs(g(i, j ) - g(i, j )) 

Where, 8{i, j) is the absolute value of difference between 
intensity g(i, j ) and its local mean value g(i, j) . 

F(i, j) is non white pixels set. It is defined as: 

Me j ) = 1 - Me j) 
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A. Neutrosophic Image Entropy: 

Neutrosophic image entropy [44] is defined as the 
summation of the entropies of three sets T,I,atldF, 
which is employed to evaluate the distribution of the 
elements in the neutrosophic domain: 



En Ns = En T + En I + En F 

maxjr} 

En T = - ^ p T (i)ln p T (/) 

i = min {T} 




Fig. 1 : Optimization image clustering flowchart. 



max {/ } 

En , = - E P/ (0 ln /^/ 0) 

/=min {/ } 
maxjF} 

En F =- YjP F (i)]np F (i) 

/=min {Z 7 } 

where En T , Ellj , and En F are the entropies of the sets T 

, I and F , respectively. p T (/), p t (i), and p F {l) are the 
probabilities of elements in T , I and F, respectively, whose 
values equal to i . 



IV. NEUTROSOPHIC IMAGE FEATURE 
EXTRACTION: 

Image feature Extraction is the first step to image retrieval 
system. Neutrosophic image (N s Im) is divided into three 

matrices are represented as images called object, edge and 
background. Each image is consisting of matrix 
representing the probability white pixel values for object 
component and probability of non white pixel values for 
background component while the intermediate matrix 
expresses the probability of the boundary between the 
white and non-white pixels. The combinations of pixel 
brightness value in(A^ s Im)components are calculated by 

using a widely method namely Gray Level Co-occurrence 
Matrix (GLCM) []. The spatially related in various 
directions with reference to distance and angular 
relationships for co-occurring pairs of pixels is one of the 
most important advantages for GLCM calculations. 
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The feature extraction lor(/V s Im) components by GLCM 

is based on pixel and its next neighbor pixel. The Contrast, 
Energy, Homogeneity and Correlation are the parameters 
of GLCM which calculated by: 

m n 

Contrast = X X (i - jf I comp (i, j ), where 

i = 1 7=1 

0 <i<m , 0 < j < n 
( m , n) is image dim entions 



Energy = comp {'h j)J 

i= 1 7=1 

m n J (/ /) 

tj v, V’V’ compV'JJ 

Homogeneity = p r 

i=i 7=1 1 + _ 7 1 

„ , (j,;)-{ A X77 } 

Correlation = > > — , where 

/=, /=, 0- A . X <T v 

fi x ,/u y are mean of probability matrix. 

cr x , cr v are 5 tan Jar d deviations of probability matrix. 

V. NEUTROSOPHIC IMAGE CLUSTERING: 
Image clustering can classify similar images into the same 

group. Let image data set be Im = \l‘ comp ,i = 1,. 2,...,n}, 

and /’ be an image in a d-dimensional space. Image 
clustering problem is to find image 
category Cl- = iCZ- , Cl- ,...,Cl ■ [, which satisfies: 

© J im t im^ ’ im 2 5 5 im n ) ’ 



Im = 0CZ„ 



1=1 

c/^. *0 fori = 1 , 2 ,..., m 

Cl imi n Cl [mj = «D /or i, j = 1,2,..., m, i * j 

Among clustering methods, the fuzzy c-means algorithm is 
widely used. An objective function for a clustering method 
is important to define. The objective function of fuzzy c- 
means is defined as: 



F* 



A(r,£)=Z2>*)" 

*= 1 k = 1 

Where Ifl is constant, and m > 1 . Cluster / is expressed 
as e i {i = 1,2 ,...,c). The membership between sample 
and cluster is expressed as; 



Eik 



f i = \,2,...,cf 



, where, 



k = 1,2,..., n y 
Ihk e {0,1}, V/, k\ =\,Vk 



can be computed by: 
Mik = ~ 



7=1 



( 


. 








Vk 


~ e i 






\x,. 


-e 




V 




7 


\) 



1 

m— 1 



where e t is the cluster center and can be computed by: 

ZU* T x k 



e ; = 



k=l 






-, l<t<c 



it=l 



The mean and the variance of ju ik for the cluster are 
computed as: 



Ei = 



a-= i 



>h 

the 



cr: = — 



>h 

partition 



matrix 



Where, the fuzzy 

is //, — , H i2 ,. . . , }■ 

As shown in figure 2, (W s Im) clustering method based on 
fuzzy c-means is used. 



A. 



Neutrosophic image clusters enhancement: 



is 



The indeterminacy of image pixel t cam ff,j ) i 

determined by its intensity value l{i,j). Strength of the 
correlation between neutosophic image components T and 
F with I are influenced by the distribution of the pixels and 
the entorpy of 7. 

The set C [o,i] may represent not only indeterminacy 
but also vagueness, uncertainty, imprecision, error, etc. So, 
the overlapping problem will appear within and between 
neutrosophic image cluster as shown in figure 3. 

Threshold processing will solve the overlapping problem 
within and between neutrosophic image clusters by 
determine the mysterious region between background and 

objects. In gray level domain, (A — mean ) operation for 
image Im GL is defined as: 



1 






Im GL (*\ J ) = X X ( m ’ U ) 



) n = J~ 



Where w is the size of the window, {m,n ) is the location 
of the pixel centered the window. A A — mean operation 



for P Ns , Pf is defined as: 



i=i 
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p A = 

1 Ns 



p(t{a)J{a),f{a)) 



, [T I <A 



T I >A 

; i w/ 

1 



j + " 



T A (Uj)= Z 



W-XW m=i ^ w/ n=j _ w/ 



F(A) 



[F 1 < A 
{F 1> A 



F A (i,j) 




j +w A 

Z F ( m ’ n ) 

"=j- w A 




Fig. 2: flowchart of neutrosophic image clustering. 
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Fig. 3: Overlapping within and between neutrosophic image cluster. 



f A (/, j)= i — E^iA 

v 7 M — M 

1 max min 

Where H{i, /) is the homogeneity value of T 7 at (i, j) . 
w is the local window size. 

After /L — mean operation subset T became more 
homogeneous after removing the noise. By using a simple 
threshold method can be segnent the subset T accurately. 

B. Optimization neutrosophic image clustering: 

Optimization (N s Im) Classification (ONsIC) method 

presents a new method to determine the best category to 
include the query images, where the characteristics of 
clusters are represented by neutrosophic sets. Suppose that 

a set of clusters Cl = {Cl , , Cl-, Cl m } which consists 



of m neutronsophic image clusters (NsIC) from which the 
most preferred cluster is to be selected to including the 
query image. Each NsIC is assessed on n different 

components as {/ m c ° mp , im c ° mp , . . . , im c ° mp } . The 

evaluation of the cluster Cl i with respect to the component 

of im c ° mp is a neutrosophic set. The neutrosophic 

index I - is such that the larger I ~ the higher a hesitation 

margin of NsIC Cl, with respect to the components of 

im , j""’ p whose intensity is given by 7T . The NsIC matrix is 

given in the following form: 



N S IC = 





i 1 

comp 


I 2 

comp 


I n 

comp 


Cl, 


C^ti Ai’^n) 


(^12’ A 2’ ^ 12 ) 


(C\n’ An’ C\„) 


Cl 2 


(^21’7u’ C 2 \) 


(^22’ ^22’ ^ 22 ) 


(T« ’ 1 2 n ’ C 2n ) 


ci m 


C^ral ’ 1 ml ’ ) 


(■^ra2 ’ 1 m2 ’ ^ m2 ) 


* ’ ^'mn ’ I mn ’ ^ mn ) 



Where the characteristics of NsIC (Cl, ) are given by: 



Obviously, 



0 <Klj+ K“ < 2 for 



Cl, 



where, l<i <m 



[ T n ,I n ,F n T i 2 ,I i 2 ,F i2 T jn ,l irl ,F irl 

N s IC(Cl ,) can present by another form as: 

cz,. = where 

is closed NsIC interval computed by: 



im""" 



e im . The 



to 



all Cl, £ Cl and 
the alternative 



degrees 

NsIC Cl, satisfies and does not satisfy the can be measured 
by the evaluation function ( E ). The evaluation function 
E(Cl , ) of alternative NsIC Cl, can be expressed as: 
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E(ci, ) = [ k'j , k; ] A [K ' ik , K“ k ] A . ... A [K ' ip , K ; J v [k! q , k“ q j 

= [min {K‘ , K\ k ,...,K < p }, mmfo , K“ k ..... K “ p }] v [k ' q , K“ Q } 

= [max [min {kI , K' ik K' ip }, k! q }, max [min {k“ , K* K‘ }, K“ g }] 
= [ K k’ K ci] 

where, 1 < i < m , and (a and v) denote the minimum 
and maximum operator of neutrosophic set respectively. 

The score of alternative NsIC S c (Cl i )can be evaluated 
by: 

SACl,)= Afa, -T ‘, J 



2^max| 


f Tc/. + 1 ci, ^ 


+M' 


| - min| 


|%, + 1 ci, 


fi-^+vdYi 


1 ^ j 


1 2 J, 


I ^ J’ 


l ^ JJJ 



An accuracy function A c is used to evaluate the degree of 
accuracy of neutrosophic elements as follows: 

\{Ch)=\(Ta,+ T a,) 



Lmff r « +/e 0 


1 a, A 


j + max| 


( + 1 a, j 


+ 

1 


■Hi 2 J’ 


{ ^ J, 


.1 ^ ) 


l 2 jjj 



The larger value of A c (Cl . ) represents the more degree of 
accuracy of an element Cl, in the neutrosophic set. Based 



on the score function S c and the accuracy function A the 



degree of suitability that the corresponding cluster Cl, 
satisfies the query images component can be measured by: 

WiEia , )) = (S c (. E(Cl ))) 2 ^ c/i ^ 



so 



The coefficients in W{e(CI , )) have been chosen 

that W{E{Cl i )) e [0,1]. The larger value of W(E(C1 , 
the more suitability to which the alternative NsIC 
(Cl , ) satisfies to query image components, where 



1 < i < m. 

Assume that there is a query images wants to choose an 
alternative NsIC which satisfies the components 
of Q1 . , where, 1 < j <n , each QI . have a different 

degree of components {t, 1 , f 'j. Where, 

0<f <1 , 

0</<l, 

0 < F < 1, 



and 0 <f + 1 + F < 3 

( T,.I,andF j are the degrees of membership (object), 
indeterminacy (edge) and non-membership (background) of 
the images un . e im to the vague concept 
importance of criterion respectively. 

The weight of query image components lies in closed 
interval [vt/ , w“- ] where. 





r fj+ij 


i-Fj+iE 




r f j+ ij 


i-Fj+i,' 


. 

min 






, max 






i 2 


’ 2 

y 


v 2 


9 2 

J _ 



Where, for each query images are 0 < w'j < w“ < 1 . 
Then the suitability degree of alternative NsIC (C/ ; ) 



satisfies the query images components can be measured by: 
«(C/,.) = max{(w(CZ,.)) 2 



i . 2 / 

The optimal weights value can be computed via the 
following programming: 

T - - cV 



i=l 

Subject to the conditions: 



1 -- 



(2 2^- 



w'j < wl < w" 



where, j = l,2,...,n 



VI. Accuracy assessment: 

A neutrosophic image is a combination of various 
components built on the intensity of pixels. The 
unsupervised classification applied to neutrosophic image 
sets result in clusters. The comparison between 
neutrosophic image components and clusters indicates the 
accuracy of the classification. Confusion matrices are used 
to assess this accuracy as they compare the relationship 
between the neutrosophic image components as the 
reference images and the corresponding results of the 
unsupervised classification technique. 

A resulting cluster by unsupervised classification is not 
automatically labeled nor identified as corresponding to 
specific image. So, image is investigated with respect to all 
clusters, and the cluster containing most of image 
components closest to the mean of that image is considered 
as its corresponding to a specific cluster. Based on 
confusion matrix, the accuracy is then expressed in terms 
of the kappa statistic (k) where the difference between the 
clustering accuracy and the chance agreement between the 
images and the clusters is calculated [45]. 

It results in a value between 0 and 1 for each classification, 
where 0 indicates that the clustering is no better than 
grouping the image by chance: 

k = ^i id 

N 2 -±(yry,) 

i = 1 

Where: 

c ■ is the number of clusters 1 < / < C 
N ; is the total number of images in the image database 
classified 

i : refers to the images corresponding to cluster i 

y • is the number of images in row i and column j in 
,J the confusion matrix 

y ■ is the number of images in row i in the confusion 
matrix 

y • is the number of images in column j in the 
1 confusion matrix 
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VII. EXPERIMENTAL AND RESULTS: 

The WBIIS image database [46] is used to evaluate the 
ONsIC system. WBIIS image database consists of 10,000 
generic images with variety size distributed on 10 



categories of Africa, Beaches, Building, Buses, Dinosaurs, 
Elephants, Flowers, Foods, Horses and Natural. Several 
test images were used during experimentation and result as 
shown in figure (4). 




Building 



Bus 







_n 52 








m 




E 


fit: 




Flower 




Food 








Horse 



Natural 



Fig. 4: Database categories 



The experimental results are discussed under 4 heading: 

A. Transform image to neutrosophic image. 

Figure 5 illustrates three steps to transformation from RGB 
images to neutrosophic images. First step convert the RGB 
image to Gray image. A gray image is a simple kind of 
image that contains one domain, and each pixel in the 



image can be represented by an integer. Second step, an 
input image is mapped to T and F by the S -function [47] 
and indeterminacy domain I by homogeneity [48], There 
are 12 features (4 parameter of GFCM x3 image 
component (T, I and F)) stored in features vectors database 
for the further processes. 




Fig. 5: Neutrosophic image transform steps. 



B. Neutrosophic image clustering. 

The major advantage of neutrosophic image clustering 
(NsIC) method based on fuzzy c-means can deal with 
indeterminate intensity regions effectively and accurately. 
Another advantage of NsIC is that it can smooth the 
complex background; therefore, it can prevent the object 



region from connecting with the false foregrounds. Besides 
the above two major advantages, NsIC finds more accurate 
object boundaries. Table 1 illustrate example of image 
clusters. 
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Table 1: Example of image clusters. 



Clusters 


#of image 


im[° mp 


im c 2 omp 


im c ° mp 


im\ omp 


T 


F 


I 


T 


F 


I 


T 


F 


I 


T 


F 


I 


Cl 

^ 1 Africa 


0.5777 


0.4223 


0.0274 


0.5503 


0.4497 


0.0209 


0.5474 


0.4526 


0.0433 


0.6357 


0.3643 


0.0280 


Cl 

Beach 


0.4135 


0.5865 


0.0183 


0.4839 


0.5161 


0.0212 


0.5977 


0.4023 


0.0131 


0.4668 


0.5332 


0.0149 


Cl 

Build 


0.4806 


0.5194 


0.0474 


0.4817 


0.5183 


0.0326 


0.4333 


0.5667 


0.0353 


0.5609 


0.4391 


0.0451 





Fig. 6: Enhanced within neutrosophic image. 



C. Neutrosophic image cluster enhanced. 

The window size (5x5) is used for computing the 
standard deviation and discontinuity of pixels intensity. 
Figure 6 illustrates the homogeneity image in domain I. 
Otsu's method [49] is utilized a simple thresholding 
method. The global t optimum threshold is finds that 
minimizes the overlap variance of the background and 
objects within cluster and the overlap variance between 
images between clusters by the following equation: 

a(t) = c ] {t)a ] {t) + C2{t)(J 2 {t) 



\K‘ ,K" = 



mil/ 1 '** ' 






T, +/„ 1-F, +/,, 



k ] 



2 2 

The result is: 

[0.7,0.725] [0.57,0.63] [0.275,0.3] [0.3, 0.6] 

" [0.6,0.675] [0.825,0.925] [0.25,0.75] [0.35,0.75]_ 

Scoring and accuracy evaluation matrices S , A c are 
computed by: 

s, = 2 x \4 h - < ] , 4 = 0.5 x [4, t + 4 h \ 

The result is: 



(f)is 


the sum of variances of the two clusters as 


S c = 


'0.05 


0.12 


0.05 0.6' 


, 4 = 


'0.7125 


0.6 


0.2875 


0.45" 






0.15 


0.2 


OO ' 
O 




0.6375 


0.875 


0.5 


0.55_ 



and probability of class i , i=l,2 respectively. Threshold t 
that results in the minimization of 0 (t) separates the two 
classes as the foreground and background, respectively. 

D. Optimization neutrosophic image clustering. 
ONsIC technique is applied on neutrosophic set of image 

clusters (NsIC) {(7/^ 7 ,/ = 1,2,..., mj with 

evaluations as neutrosophic components. For example, 
there are two clusters each cluster including four images 
components (T, I, F) are represent as: 

/' I 2 1 3 I 4 

comp comp comp comp 

N s C = Cl l [0.8,0.6,015] [0.68,0.46,0.2] [0.45,0.1,0.5] [0.5,0.4,0.8] 

CU [0.4,0.8,0.45] [0.75,0.9,0.05] [1,0.5, l] [o.5, 0.6. 0.9] 

Evaluation of neutrosophic component for NsC 

(e(MsC ij )) = \k\j , ] is computed by from: 



Weights of NsC W{NsCl ij ) are computed by: 

w[Nsci ij )= (s c {mc, ,)} - 1 ~ 4 ( iv ' c ^ 



The result is: 

w{n s ci,j)= 



-0.141 -0.186 -0.354 0.085 
-0.159 -0.022 0.75 0.415 

The coefficients of the linear programming problem are 
computed by sum each column in W^NsC^ ) : 

XwkaJ=[-0.3 -0.208 0.396 0.5] 

Four neutrosophic query images 

1 QHfj.ij.Fj) ’ J = are given as: 
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comp 

n Q i 



comp 

n Qi 



QI (t,i,F) |^o.25,0.3,0.25 ’ 0.35,0.6,0.41 ’ 0.32,0.55,0.67 ’ 0.64,0.98,0.57 I 

For simplifying computation, the neutrosophic set may be 
written as: 









im c n omp 

til 3 


im c ° mp 

til 4 


T 


0.25 


0.35 


0.32 


0.64 


I 


0.3 


0.6 


0.55 


0.98 


F 


0.25 


0.41 


0.67 


0.57 



Weights of query image are computed by formula: 



T, + / 



l — F.+I, 



ff 



T,+I , 



l — F.+I, 



JJ 



The result is: 

k.wjJ =110.275,0.525] [0.475,0.595] [0.435,0.44] [0.705,0.8 1] 
Thus the linear programming now can be set as: 

max Z[Z w ( A “)j 

Subject to: 



Wj < w J < w'j 



The result is: 

Maximize : 
subject to : 



-0.3w, -0.208 w 2 +0.396w 3 +0.5w, 
0.275 < w, < 0.525 
0.475 < w 2 < 0.595 
0.435 <w,< 0.44 
0.705 < w 4 <0.81 

The linear programming can easily solve to get wl : 



wi = 0.275 
wl = 0.475 
wl = 0.44 
wl =0.81 

Then the degree of suitability that the corresponding cluster 
Cl i satisfies the query images can be measured by the 
following function: 

(4k -K f t 2 — 



«(«,)= z 

j = 1 



The result is: 

R(c/,) = ^0.05f 

+ (Vo5) 2 - 

R(C/ 2 ) = f(0.15) 2 - 



1-0.7125 



1-0.5 

2 



0.44 
1-0.6375 



X 0.275 + 1 (0.12)' 

1-0.45 



1 - 0.6 

2 



X w* 



< 0.475 



(0.6) 2 



x 0.275+ (0.2) 2 - 



x 0.81 = -0.21380375 
1-0.875 



+ [(l) 2 -[^)]x0.44 + ^(0.8) 2 - 



1-0.55 

2 



x 0.475 



x 0.8 1 = 0.6 11 80625 



Therefore we can see that the alternative cluster Cl 2 is the 
best choice. And the optimal ranking order of the 
alternatives is given by Cl 2 {Cl x . 



VIII. Performance results 

Table 2 is show the confusion matrix for NsIC technique, 
the ten clusters found that represent a good approximation 
of the ten categories of database image. 

Table 2: Confusion matrix for NsIC. 



Categories 


Clusters 


a 


C/ 2 


c/ 3 


CI t 


Cl s 


ci t 


Cl. 


c/, 


c/, 


CL 






0 


5 


0 


1 


0 


3 


0 


2 


0 




0 


98 


1 


0 


0 


0 


0 


0 


0 


i 


Cat^ 


2 


0 


21 


0 


1 


0 


1 


0 


5 


0 




0 


0 


5 


87 


0 


6 


0 


1 


1 


0 




0 


6 


0 


1 


92 


0 


1 


0 


0 


0 




6 


2 


0 


0 


0 


88 


2 


1 


1 


0 


Cal r,~ 


0 


0 


0 


2 


0 


3 


95 


0 


0 


0 


Ca! r,~ 


0 


0 


0 


1 


2 


0 


0 


21 


0 


0 


CtA 


0 


3 


0 


0 


0 


0 


1 


0 


93 


3 


Car. vi w 


3 


0 


1 


0 


1 


0 


0 


0 


5 


90 



The experiment results have proved beyond any doubt that 
NsIcC technique can be of high-accuracy techniques. 



Where, the accuracy of similar images was gathered in 
different clusters exceeded 90%, as shown in figure 7. 
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Fig. 7. The accuracy assessment (k) versus different clusters. 




Fig. 8. The recall and precision measures versus different clusters. 



The retrieval performance is defined by the precision and 
recall of the retrieved images [50]. As shown in figure 8 the 
beach category has been identified with high recall 98%. 
That is mean in this category the high number of similar 
images are gathered with very few errors. Whereas the low 
recall 80.8318% has been identified for house category. 
That is meaning there are misclassified for some images 
in this category as a result of the weakness of I set that is 
working to determine the threshold used in image 
enhanced. The precision value ranging from 86.9195% to 
97.9798% where the highest value belonging to the cluster 
8 and lowest for the cluster 9. This means there is a cluster 
8 has more than an image that dos not belong to texture 
category and vice versa in cluster 5. 

The experiments were carried out on a large-scale over the 
WBIIS image database to measure the effectiveness of 
NSIC. The results focused on recall and precision measures 
and accuracy assessment as a criterion for comparison 
between NsIC technique and other technique namely G. 
Ciocca et al (2013). 

G. Ciocca et al [51] proposed a new technique for 
unsupervised image classification based on purely semantic 
representation not on statistical method. This technique is 



composed of two parts namely supervised feature [52] and 
unsupervised classification. The idea of supervised features 
part is based on the hypothesis in the semantic realm 
regarding categories. That an image d belongs to category 
C i with certain probability will not determine whether the 

same image belongs to a category C j , but it will modify the 

probability that it be so. There are supervised feature 
descriptions are used to classify image regions into 
semantic concept classes such as Classemes [53], 
Prosemantic features [54], Object bank [55], CCA-based 
features [56], and Primitive features [57]. The second part 
unsupervised classification algorithm such as k-means 
based on the output of a limited number of classifiers can 
be embedded into a feature space by a clustering algorithm. 
G. Ciocca et al is used prosemantic features as supervised 
feature and k-means algorithm as unsupervised classifier 
and applied on image database. The results of recall, 
precision and accuracy assessment indicates that the NsIC 
technique achieved a high performance compared with the 
other technique as shown in figures 9, 10 and 11. 
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Fig. 9. The recall for compared two techniques. 



In general, it is obvious that the NsIC technique has 
achieved a better precision, recall and accuracy assessment 
(k) at various numbers of clustered images than the other 
compared technique. 

IX. CONCLUSION: 

This paper presents a new technique to unsupervised 
classification for images based on neutrosophis sets and 



optimization linear programming. Neutrosophic theory was 
used to transform the gray image to neutrosophic image 
components (O, E, B). Indeterminacy set (E) was worked 
on determine the objects boundaries with high precision. 
Determining the boundaries of objects accurately blunted 
the effect of the overlapping problem within the cluster. 




Fig. 10. The precision for compared two techniques. 




Fig. 1 1 . The k statistic for compared two techniques. 



Neutrosophic image clustering method based on fuzzy c- 
means is used. Neutrosophic image clustering has been 
enhanced by using the /.-mean operation which helped on 



solve the overlapping problem between clusters. 
Optimization neutrosophic image clustering is achieved by 
using the weight coefficient between image clusters and 
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images category as an object function in linear 
programming problem. Whereas, the constraints of linear 
programming problem are the weight limits for query 
images. 

Practical results conducted on neutrosophic image 

clustering technique has proved its efficiency where it was 
to obtain the high performance rate in the accuracy of the 
resulting clusters as well as high values of recall and 
precision measures. 
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